II 



-Loan copy: return to 

AFWL TECHNICAL LIBRARY 
KIRTLAND AFB, N. M. 


KALMAN FILTER ESTIMATION 
OF HUMAN PILOT-MODEL PARAMETERS 


James R. Schiess and Vincent R. Roland 

Langley Research Center 
Hampton, Va, 23665 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • NOVEMBER 1975 


TECH LIBRARY KAFB, 



1 Report No 

NASA TN D-8024 

4. Title and Subtitle 


2. Government Accession No. 


KALMAN FILTER ESTIMATION OF HUMAN 
PILOT-MODEL PARAMETERS 


7. Author(s) 

James R. Schiess and Vincent R. Roland 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 


12 Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

15. Supplementary Notes 


tech library kafb, nm 



3. Recipient's Catalog No. 

5. Report Date 

November 1975 

6. Performing Organization Code 

8. Performing Organization Report No. 

L-10041 

10. Work Unit No. 

506-25-99-01 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

14. Sponsoring Agency Code 


16 Abstract 


The parameters of a human pilot-model transfer function are estimated by applying 
the extended Kalman filter to the corresponding retarded differential-difference equations in 
the time domain. Use of computer-generated data indicates that most of the parameters, 
including the implicit time delay, may be reasonably estimated in this way. When applied 
to two sets of experimental data obtained from a closed-loop tracking task performed by 
a human, the Kalman filter generated diverging residuals for one of the measurement types, 
apparently because of model assumption errors. Application of a modified adaptive technique 
was found to overcome the divergence and to produce reasonable estimates of most of the 
parameters. 


17. Key Words (Suggested by Author (s) 1 


18. Distribution Statement 


Kalman filter 
Pilot model 
Adaptive technique 


Unclassified — Unlimited 


Subject Category 64 


19. Security Classif. (of this report) 

1 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

42 

$3.75 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


KALMAN FILTER ESTIMATION OF HUMAN 
PILOT-MODEL PARAMETERS 

James R. Schiess and Vincent R. Roland 
Langley Research Center 

SUMMARY 

The parameters of a human pilot-model transfer function are estimated by applying the 
extended Kalman filter to the corresponding retarded differential-difference equations in the 
time domain. Use of computer-generated data indicates that most of the parameters, includ- 
ing the implicit time delay, may be reasonably estimated in this way. When applied to two 
sets of experimental data obtained from a closed-loop tracking task performed by a human, 
the Kalman filter generated diverging residuals for one of the measurement types, apparently 
because of model assumption errors. Application of a modified adaptive technique was found 
to overcome the divergence and to produce reasonable estimates of most of the parameters. 

INTRODUCTION 

The mathematical modeling of human dynamics in specific tasks is of considerable 
interest to control engineers. For example, a model which accurately predicts a pilot’s 
response in the pitch control of a spacecraft is beneficial in the analysis of the handling 

qualities of the vehicle during the development stages of the vehicle. Studies of the pilot 

models in general have been reported in the literature (refs. 1 to 4). 

In order to use any of these models, it is necessary to know various parameters which 
occur in the model. Investigations of applying model-matching techniques to determine the 
best parameters from pilot-response data have been reported previously (ref. 5). Since the 
human pilot-model equations classify as retarded differential-difference equations in the time 
domain, the application of linear estimators to this type of equation by Koivo and Stoller 

(refs. 6 and 7) and by Kwakernaak (ref. 8) is also of interest. The work by Kwakernaak 

(ref. 8) on the linear problem concerns an estimate of the state only and no parameters. 
Although an exact solution for the optimal filter is obtained, this approach presents practical 
computational difficulties. The findings of Koivo (ref. 6) present an estimator (not neces- 
sarily optimal) both for the state and for particular parameters, but does not include the 
time delay. 

The present paper investigates an application of a form of the extended Kalman filter 
to the time-domain representation of a particular pilot model in a single-control compensatory 



tracking task. The time delay is included as one of the parameters to be determined. The 
time delay is retained as an implicit parameter rather than as an explicit parameter (obtained 
if a Fade approximation is applied) for exactness. When applied to the experimental pilot- 
response data, the initial results of the Kalman filter indicate difficulties in obtaining satis- 
factory estimates. An analysis of the mathematical model and techniques for improving the 
estimates are presented. 

SYMBOLS 

c(t) plant output, V 

D(t) input-disturbing function, V 

E{| expected value of | | 

e(t) system error, V 

f(y(t),y(t - 0),) 13 X 1 column vector representing right-hand side of differential- 


’ D(t - 9),t) 

) difference equations 


h('k-h) 

measurement coefficient matrix at tj.^j 


I 

identity matrix 


K 

plant static gain, sec' 


*^1 

pilot static gain 


K(tk+l) 

Kalman filter gain at tj^.|.j 


P(tk+l|^k) 

state covariance matrix at tj ^.(.2 before processing measurement at 

h+l 

P(tk+1 jtk+l) 

state covariance matrix at t]^.j.j after processing measurement 

at tj,+i 

Q(^k+l) 

process noise covariance matrix 


q 

integer defined by q = 0/At 
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q(k) 

qOc + 1) 

^(^k+l) 

Rp(s) 

S 

t 

At 

u(t) 

w(t) 

Wp(t) 

Yp(s) 

y(t) 

V4(t) 

" l>«2 

r(k + 1) 

S(t) 

T?(t) 

d 


variance computed in adaptive procedure 
estimated process noise variance at tj^+j 
measurement covariance matrix at tj^_j_j 
second-order transfer function in pilot remnant 
Laplace operator, sec"^ 
time, sec 

k plus first time t, sec 
time increment, sec 
deterministic pilot function, V 
13 X 1 process noise vector 
white noise input, V/sec^ 
deterministic pilot transfer function 
13 X 1 state column vector 

fourth component of y(t) defined in equation (6), V 

coefficients in pilot remnant transfer function 

process noise coefficient matrix at 

stick motion, V 

pilot remnant, V 

pilot time delay, sec 



pilot lag-time constants, sec 




pilot lead-time constant, sec 

0/tj^-|.j,tj^j state transition matrix from tj. to tj^^_j 


Superscripts: 


^ best estimate 

• first derivative with respect to time 

second derivative with respect to time 

/ 

T matrix transpose 

-1 matrix inverse 

Subscript: 

k denotes evaluation at time tj^ 

PROBLEM DESCRIPTION 

This analysis considers the problem of identifying the parameters in an assumed input- 
output relationship describing the human pilot from data measured in a tracking task in which 
the pilot provides feedback control. Because of process noise which is represented by the 
pilot remnant and measurement noise in the data accumulation, the estimation process con- 
siders the problem of minimizing the effect of the noise. A well-established method of 
reducing the error in such an estimation process is the Kalman filter. This estimator con- 
siders both process and measurement noise and provides the best estimation of state and 
parameter vectors in the sense that the expected value of the sum of the squares of the 
error in the measurement is minimized. In addition, this filter provides a covariance matrix 
which indicates the quality of the estimates. 

A block diagram of the pilot airplane model used in this analysis is given in figure 1. 
The model consists of a closed-loop single-control compensatory tracking task. The dynamics 
of the controlled elements are of the acceleration type (k/s^|. Such a type is found in the 
attitude control of a space vehicle by control jets. The dashed lines enclose the portion of 
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the model representing the pilot. The transfer function representing the deterministic portion 
of the pilot model Yp(s) is given (ref. 4) as 


Yp(s) = 


u(s) ^ 
e(s) 


(rjS + l)(,^s f 1) 


( 1 ) 


The parameter is a pilot static gain and 6 is an effective time delay made up of 

transport delays and high-frequency neuromuscular lags. The constant is a time con- 

stant associated with the low-frequency characteristics of the neuromuscular system (ref. 4).i 
The parameters and are pilot-adjusted parameters. 


The pilot-remnant function used in this analysis is a second-order noise filter which 
converts the white noise input Wp(t) to colored noise input Tj(t). The form of this 
filter is 


Rp(s)> , ' — (2) 

s + 

The remnant 77 ( 1 ) represents noise produced by the pilot. Such noise is caused by high- 
frequency neuromuscular lag and lead terms and low-frequency neuromuscular lag terms. 

For the current analysis the measurement types used are the plant output c(t) and 
the stick motion 6(t). The latter is the sum of the deterministic pilot function and the 
pilot remnant. 

The functions of time which represent the aircraft dynamics, the deterministic pilot- 
transfer function, and the pilot-remnant function, respectively, are given by the differential- 
difference equations 

c(t) = K 5(t) (3) 

u(t) + u(t) = K 

■j7(t) + ^(t) + 7?(t) = Wp(t) (5) 

These equations indicate two possible difficulties in applying the Kalman filter; (1) because 
of the way in which 6 occurs, they represent a hereditary system (i.e., one represented 
by differential-difference equations); and (2) the equations are nonlinear in the combination 
of the variables and parameters. The problem of nonlinearity can be overcome by applying 


T r 
1 2 


u(t) + 


(h 


+ r 




V(t 


0) + e(t - 0) 


(4) 
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the extended Kalman filter to the nonlinear equations. The problem of applying the Kalman 
filter to a hereditary system is considered in the section on “Filter Equations.” 

For convenience in the application of the Kalman filter, equations (3) to (5) are trans- 
formed to a state space representation by using the following definitions for the time- 
dependent variables y^^ft) to yg(t)- 


yj(t) = c(t) 


y^Ct) = c(t) 


y^ft) = u(t) 




K 


( 6 ) 


y^(t) = u(t) + u(t) - — r^e(t - 6) 


y^(t) = T?(t) 


y,(t) = 7?(t) 

O 




The variable is defined in this way to simplify the final state equations. With these 

definitions and with the use of the fact that e(t) = D(t) - c(t), the state space form of 
the pilot-model equations is given by equations (7) to (12) 


Yj(t) = 

y^(t) = 
^ 3 «> = 
y^(t) = 

y,(.) = 

y^(t) = 

o 


y,(t) 

K[y3(.) + y^Ct)] 

J-y(t)-^y(t)y^TD(t-«-y(t-e)1 
2 2 1 2 L -J 

K / ’■ \ 

■ -r y/t) + - d) - y (t - 0)] 

1 ^ 1\ 1/L -» 

yit) 

O 

-a y,(t)-a y (t) + w (t) 

1 O Z J ^ 


(7) 

( 8 ) 

(9) 


( 10 ) 


( 11 ) 

( 12 ) 
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It should be observed at this point that equations (7) to (12) are linear in the y’s. 

Thus, if y(t) is defined as the six-dimensional column vector whose components are 

to yg(t), A is the 6X6 constant matrix containing the coefficients of the right- 


yj(t) 

hand sides of these equations, 


b is the six-dimensional column vector containing the coef- 


ficients of D(t - 0) - y^Ct - 0), and w(t) is the vector of white noise input. Equa- 
tions (7) to (12) can then be written in vector form as 


y(t) = Ay(t) + b|D(t - 0) - y^(t - 0)J + w(t) (13) 

Equation (13) is in the form of a single-input, time-invariant linear control system to which 
the Kalman filter is immediately applicable, (See, for example, ref. 9.) 

In addition to the six variables of equations (7) to (12), it is also of interest here 
to estimate seven of the parameters in these equations. The estimation may be accomplished 
by augmenting the state equations with the system 


yi(t) = 0 


(i - 7, 8, . . ., 13) (14) 


where y^(t) through defined as the parameters K^, r^, 0, 

and o:^, respectively. With y(t) extended to a 13-dimensional vector with the newly 
defined components, the right sides of equations (7) to (12) and (14) are now nonlinear 
functions of the vector y(t). The general form of these equations is given by 

y(t) = f[y(t), y(t - 0), D(t - 0), tj + w(t) (15) 

where 1 [y(0’ y(l ' 0)^ ^^(1 ‘ 0)^ l] ^ 13-dimensional column vector defined by the 
right sides of equations (7) to (12) and (14), and w(t) is a 13-dimensional column vector 
which is all zero except for the sixth component which is Wp(t). The time delay 0 is 
left explicit in equation (15) for later reference. 

In the vector notation, the measurements are given by 


m(t) = H(t) y(t) + V(t) 

where m(t) is a two-dimensional column vector. 


fo 0 1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

o’ 

O 

O 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


(16) 


7 


II 



where V^(t) and ^2^^^ represent white noise and are assumed to be uncorrelated with 
each other and with the process noise Wp(t). In equation (16), tht first measurement 
is 5(t), the second is c(t). 

In subsequent sections of this report the stability of the closed-loop system is discussed. 
Because the time delay occurs implicitly in equations (7) to (12), the effect of particular 
values of this parameter on the system stability cannot be readily discerned. Instead, if the 
first order Fade approximation to die exponential is used in equation (1), and if equa- 
tions (2) and (3) and the new equation (1) (resulting from use of first-order Fade approxi- 
mation) are transformed to the time domain, the resulting equations are of the form 

y(t) = A' y(t) + b' u(t) + w(t) 

where 6 occurs explicitly in the new A' and b The approximate asymptotic stability 
of the system can then be determined by examining the eigenvalues of A'. Thus, in this 
report, the system is considered stable if all the eigenvalues of A' have negative real parts. 

FILTER EQUATIONS 

The filter equations used for this study are the discrete extended Kalman filter equa- 
tions as given by Jazwinski (ref. 9). The use of these equations overcomes the difficulty 
which arises from the fact that the state equations are nonlinear when the parameters are 
treated as state variables. It should be pointed out that the extended Kalman filter 
actually is a nonlinear filter; for example, the Kalman filter is identical to the nonlinear 
discrete invariant embedding algorithm of Sage and Melsa (ref, 10) when the measurement 
equation is linear as in the case under consideration here. 

In general, with any discrete filter the state is predicted to time tj^_j_| through the 
application of the state equations to the best estimate of the state at the previous time tj^. 
Since the pilot-model equations contain the state at tj^ - 0, an appropriate value of the 
state at that time must be chosen. The choice is made by assuming that the time incre- 
ment At between measurements is constant and that the delay 0 is an integral multiple 
of this increment. Thus, let q be the integer defined by 
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Then the state is evaluated at time \ ^ t>y using the best estimate of the state for the 

time tj^_q ~ ^ Since the state at is known at tj^, the delayed state can 

be treated simply as a known forcing function in the state equations. In a notation similar 
to that of Jazwinski (ref. 9) the extended Kalman filter is given in equations (18) to (22) 
with q defined by equation (17) 

y(tk+i|tk) = y(tk|tk) +/*'''^*l|'y(tkhk) y(tk-q|tk-q) ^(tk-q)’ 

^k ^ 

^(^k+l hk) “ *^^k+l’^k^ ^^k|^k) *^^^k+l>^k^ Q^k+1^ 

K(tk+l) = p(k+l hk) H'T(tk+l) [“(k+l) p(tk+l hk) H'T(tk+i) - R(tk+i)j 

y(k+l hk+l) = y(k+l hk) ^(k+l) [lll(k+l) ' P^(k+l) y(^k+l hk)j 

X 

P('k+1 hk+l) = [■ - K(tk+l) H(tk+l)j P(tk+I |lk) [> - K(tk+l) H('k+l)] 

+ K(tk+i) R(tk+|) KT(t|^^,) 

In these equations, y(^k+lhk) estimate of the state at tj^_|_j based on 

measurements through time tj^; is the covariance matrix of y^^k+l hk) 

K(tk+j) is the gain matrix. The matrix 0(lk+l’^k) transition matrix which is 

obtained by integrating with respect to time the state equations (eq. (15)) differentiated with 
respect to the state variables. The matrices Q|tj^^_jj and K(tj^+l) are the diagonal covari- 
ance matrices of the process and measurement noise, respectively. 

The effectiveness of the filter equations in estimating the state variables can be deter- 
mined by examining the diagonal elements of the covariance matrix P(tkl^k)’ these 

elements are the variances of the estimated state variables. Further, the sum of the squares 
of the measurement residuals may be treated as a cost function in order to measure the 
overall effectiveness of the filtering process. 

RESULTS FROM COMPUTER-GENERATED DATA 

In order to observe the behavior of the pilot-model state equations and to uncover 
difficulties associated with applying the extended Kalman filter to these state equations, the 
Kalman filter was first applied to data generated with the differential-difference equations. 
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(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 
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All of the computations required for this study, including the generation of the simulated 
data and the filtering of both the simulated and experimental data, were performed on the 
Control Data series 6000 series computers at Langley Research Center. 

In generating the data, reasonable initial conditions for the state variables and parameters 

were chosen; the state equations were numerically integrated to obtain a “true” time history 

of the variables. The “true” variables were then used to calculate a series of measurements 

2 

containing an additive random noise of zero mean and a variance of 0.00092 volt . For 

-9 

the simulated data presented here, the static gain K was 10 seconds , and the variance 
of the process noise was 59 600 volts^/seconds'^. The value of 59 600 volts^/ seconds^ for 
the process noise was obtained experimentally so that the ratio of the mean-squared pilot 
remnant to the mean-squared stick motion is about 0.5, as is indicated by current pilot-model 
theory (ref. 5). The parameter values along with the initial values of the variables are given 
in the first column of table I; for these values, the system is stable. (The first-order Fade 
approximation is stable.) 

The disturbing function D(t) was chosen to be 0.25 sin t. The state integrated to 
the final time of 10.25 seconds is found in the first column of table 11. Using these data, 
measurements containing the additive noise were generated at the rate of 100 measurements 
per second. 

The values used to initiate the estimation procedure are given in the second column 
of table I. The parameter values are simply guesses; the variable values were mainly derived 
from the data according to the following scheme: the second measurement (the plant out- 

put c(t)) gives a direct measure of the first variable; since the first measurement (the stick 
motion 6(t)) is the sum of the variables u(t) and T7(t), the variable u(t) was arbitrarily 
chosen to be zero and 7?(t) evaluated with the measurement value; the second variable 
c(t) is evaluated by computing the first-order divided difference of the plant output c(t) 
based on two consecutive values of c(t); the fourth variable was arbitrarily set to zero and 
the sixth evaluated by using the first-order divided difference of the measured stick motion 
based on consecutive values of 5(t). 

For the parameter variances, each standard deviation was set equal to 10 percent of 
the corresponding parameter value. Since the plant output is both a state variable and a 
measurement, it is reasonable to choose the initial variance of this variable to be equal to 
the variance of the corresponding measurement (0.00092 volt ). Since the stick motion 
is the sum of the deterministic pilot output and the pilot remnant, it is true that for the 
variances 


Ej[u(t) + 77(t)]2j = E{[6(t)]2j = 0.00092 volt^ 


(23) 
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However, expanding the left expectation yields 


Ej[u(t) + Tj(t)]2j = eju^Ct) + 2u(t) r?(t) + 772 (t)| 


= Elii2(t) + 2E u(t) T?(t) +EU2(t)| 


where in the last step the linearity of the expectation has been used. If, to simplify the 

analysis, it is assumed the two variables are uncorrelated, then the second term on the 

right in equation (24) is zero. Finally, by assuming the two variable variances are equal 

and by using equation (23), the result is that the initial variances for u(t) and rj(t) 

2 

are 0.00046 volt . The remaining variances were arbitrarily chosen to be 0.1. 

The results given in table II as case 1 indicate that c(t) and c(t) have been well 

estimated but that the remaining variables are in error by at least 8 percent. The estimates 

of u(t) and ? 7 (t) illustrate that although the error in the stick motion (which is the 

sum of these two variables) is small, the errors in the individual estimates may be large 

offsetting values. Further, an examination of the correlation coefficients, which can be 

derived by normalizing the covariance matrix, indicated that the correlation coefficient 

for u(t) and 7^(0 is -0.991. On the other hand, the error in Y^(0 is quite large; 

this large error presumably results from defining in terms of u(t) (see eqs. (6)) 

which contains a 46-percent error. Finally, the variable ? 7 (t) is poorly estimated since 

in the state equations this variable is a function of r/(t) and of the two poorly estimated 

parameters a and a . 

1 2 

An examination of the state equations (eqs. (7) to (12)) shows potential difficulties 
in estimating the parameters. These difficulties are verified by the case 1 results. First, 
the pilot static gain Kj always occurs in combination with and Although 

these latter two parameters occur separately elsewhere, the manner in which Kj occurs 
suggests difficulty in separating the three parameters. The case 1 results show that the 
estimates of these three parameters are less accurate than the initial guesses; the results also 
indicate that separation of these parameters may be difficult. Also, the delayed term coef- 
*^ 1^3 

ficients and 

^T2 

the poor estimates obtained for the parameters. 

In the case of time delay Q the estimate has changed only slightly, possibly 
because of the poor estimates obtained for the coefficients of the delayed terms in equa- 
tions (9) and (10). There is also only a slight change in The poor estimate obtained 

for Q!j is evidently because of the poor values obtained for T 7 (t) and ri(t); all these 
values occur in equation (12). 


are less accurate than the initial guess as a result of 
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In table III the variable and parameter standard deviations and the root-mean- 
squared (rms) measurement residuals are presented. For case 1 the standard deviations for 

and have changed only slightly from the initial values (0.31623). This fact, 

along with the small changes in the estimates of ocj and suggests that these parameters 

are not “observable” in this particular problem; that is, it is not possible to determine the 
values of these parameters by observing the behavior of the measurements. However, actual 
proof of observability in nonlinear problems is difficult to obtain. Instead, a separate case, 
identical to case 1 except that and were not estimated, was attempted. The 

resulting estimates for all the other variables and parameters were not significantly different 
(less than 0.8 percent) from those given for case 1 in table II; thus, retaining aj and «2 
as constant unestimated but erroneous values does not significantly bias the other estimates. 
Therefore, and 0^2 need not be estimated; however, in all the cases studied these 

parameters were estimated in order to retain any influence they may have exerted on the 
other estimates. 

One other notable item given in table III for case 1 is the rms residual for the stick 
motion 5(t). The rms residual should be approximately equal to the standard deviation of 
the noise added to this measurement (0.03033 volt). The smaller value obtained by the 
Kalman filter is undoubtedly a result of the fact that this measurement is the sum of the 
variables u(t) and ? 7 (t). Since these variables have nearly equal variances and a large 
negative correlation, the correlation tends to produce the small measurement variance. This 
tendency is clearly seen in the following relationship where Var(x) - Variance of x (any 
variable) and p is the correlation coefficient: 

1/2 

Var[u(t) + T 7 (t)] = Var[u(t)] + Var[77(t)] + 2p{Var[u(t)] Var[77(t)]} 

Another possible reason for the poor parameter estimates is the standard deviation 
of r?(t) which is much larger than the standard deviations of the parameters; this 
relationship is seen in table III where the standard deviations at t = 10.25 are given. 

The large standard deviation for r?(t) is a direct result of the large process-noise standard 
deviation. Since the variance of T](t) is so large, this variable seems to absorb too much 
of the stick-motion residual. In order to test this possibility the initial parameter standard 
deviations were increased to 100 percent of the parameter values in case 2. 

The results of case 2 in table II indicate some improvement in the estimates of 
Kp T 2 , Ty and \ however, the estimates of Tp 0, and are less accurate. 

Further, the coefficients of the delayed terms are not significantly better, evidently because 
of the degraded time-delay estimate. Therefore, increasing the parameter variances generally 
did not improve the estimates. 
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Additional experiments in choosing the variances were attempted. None of the experi- 
ments resulted in generally good estimates. Although using a smaller process-noise variance 
produces a smaller mean-squared pilot-remnant to mean-squared stick-motion ratio, the 
remaining simulated cases discussed here use data generated with a smaller process-noise 
variance. In this manner the effect of a smaller variance on the estimates was investigated. 
Except for this smaller variance, the initial parameter and variable values are identical to 
those used before. For case 3 the process-noise variance was chosen to be 576 volts^/sec- 
onds^; compared to the previous data, the ratio of the mean-squared pilot remnant to the 
mean-squared stick motion for the new data is only 0.04. 

The results given in tables II and III as case 3 were produced by using 10 percent 
of the parameter values for the initial standard deviation, extending the time interval for 
processing the measurements to 15 seconds, and iterating the filter process three times over 
this interval. Of the various attempts at filtering this data, case 3 gave the best estimates 
of the parameters. The variable estimates are not given in table II since a change in the 
process-noise variance significantly changed the time histories of the variables. However, 
the estimates of all the variables are better than in case 1 or 2; for example, the error 
in u(t) is only 10 percent. The estimates of K^, and are better than 

the previous two cases. Although the estimate of is poor, the improvement in the 

other three deterministic pilot parameters resulted in much better estimates of the delayed 
term coefficients. The degraded estimate of T 2 emphasizes the apparent difficulty in 
separating these four parameters. 

The small changes in d and 0:2 suggest that these variables are not observable for 
some reason. An “observable” parameter here means that a change in the parameter results in 
an observable change in one or more of the measured quantities. The use of large parameter 
variances in a case not presented here did produce a significant ciiange in the estimated time 
delay thus indicating that 6 is to some extent observable. However, for the same case, 
the estimate of 0^2 was not significantly different from the value given in case 3. More 
significant, perhaps, was the fact that the resulting parameter estimates yielded an unstable 
system. As a final test of the observability of 0, a new set of measurements was generated; 
this set differed from the data for case 3 in that only the input disturbance D(t) was 
changed. For the new data, D(t) was composed of two sine functions having distinct 
amplitudes and frequencies. The purpose of this test was to determine if the observability 
of 9 depends on the characteristics of the disturbing function. The time delay estimated 
from the new data (0.0402) differed only slightly from the case 3 estimate. Thus, the 
ability to estimate 9 depends on other factors. The observability of 0^2 questionable. 

Figures 2 and 3 present the time histories of the stick-motion and plant-output 
residuals for case 3. Although all the parameters are not correctly estimated, these figures 
indicate that the measurement residuals are generally random as is desired. In the case of 
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the plant output, the residuals in figure 3 are of the same magnitude as the noise; this 
finding is further exemplified by the rms residual given in table III. In addition, figure 2 
and the rms stick-motion residual in table III illustrate the fact that the estimates of 
i 7 (t) and u(t) contain a significant portion of the measurement noise. 

'y 

For the final case presented, the process-noise variance was reduced to 25 volts^/sec- 
onds^. The object was to process data with process noise comparable to that of the experi- 
mental data, even though the ratio of the average pilot remnant to the average stick motion 
is less than 0.01. Other than this smaller variance, all the parameter and variable values and 
the measurement-noise variance are the same as in the previous cases. The disturbing func- 
tion is the same function used in case 3. The estimates given in table II as case 4 are 
the result of iterating the filtering process three times. Since this estimation problem is 
nonlinear, it is reasonable to iterate the estimation process in order to obtain satisfactory 
parameter estimates. In this context “iteration” means using parameter estimates and the 
corresponding variances at the final time of one pass through the filter equations as initial 
values for the next pass through the filter. 

The results of case 4 show improvement in the estimates of Tj, 0, and 

ay In particular, the much better estimate of 6 indicates that the ability to estimate 
this parameter is mainly related to the magnitude of the process noise and not to the 
disturbing function. Further, the improved estimate of a^ is also evidently a result of 
the smaller process noise, whereas the slight degradation of the 0:2 estimate suggests this 
parameter is unobservable. The unobservability is further confirmed by the fact that a^ 
has the smallest correlations of all the parameters (specifically, the correlations between the 
parameters and the three variables (c(t), u(t), and ? 7 (t)) defining the measurements). The 

degraded estimate of T 2 indicates the apparent difficulty in separating the four deterministic 
pilot parameters. Even with this poor estimate, the two resulting delayed-term coefficients 
are generally good. 

Comparison of the standard deviations in table III with the estimation errors derived 
from the estimates in table II provides a means for determining the ability of the standard 
deviations to predict the actual variable and parameter errors. In almost every instance the 
estimate errors are less than three standard deviations; in case 4, which is of major interest 
here, over half of the estimate errors are less than one standard deviation. The one excep- 
tion to this statement is the parameter T 2 which was poorly estimated. Therefore, the 
standard deviations give a good measure of the errors in the estimates. 

Based on the results obtained by processing simulated data, certain conclusions have 
been reached. Reasonably accurate estimates of the variables c(t), u(t), and 77 ( 1 ) may 

be obtained, although separation of the stick motion 5(t) into good estimates of u(t) 
and ? 7 (t) may be difficult. If the process noise is approximately the same order of 
magnitude as the pilot-remnant parameters a^ and a 2 , then good estimates of all the 


14 



II I nil 


I III I I I 


parameters (except and possibly one of the deterministic pilot-model parameters) may 

be obtained. Further, the ability to estimate the time delay depends on the magnitude of 
the process noise and not on the disturbing function. Finally, the standard deviations pro- 
vided by the filter give a good measure of the estimate errors. 

RESULTS BASED ON EXPERIMENTAL DATA 

The experimental data used in this analysis were obtained from closed-loop compensatory 
tracking tasks conducted by Langley Research Center personnel with engineers and test pilots 
used as subjects. A description of this research has been reported by Adams (ref. 5); in the 
particular cases discussed here, an engineer was the subject. ' 

The variances of the measurements used in this experiment were obtained by setting 
the disturbance D(t) to zero, removing the pilot from the loop, and recording the measure- 
ments (which therefore should consist only of noise). The resulting value for both measure- 
ment variances was 0.00092 volt^. The variance of the process noise was obtained by a 
power spectral-density analysis of a model of the pilot remnant used by Adams (ref. 5). 

The resulting variance values were 19.8 volts^/seconds"^ and 5.9 volts“/seconds^, respectively, 

for the first and second sets of experimental data presented here. The same engineer was 
the subject for both experiments, 

In order to obtain initial conditions for the estimation procedure a modification of 
the method presented in the previous section was used. The variables c(t) and c(t) were 
evaluated from the measurements as previously described. However, the pilot remnant was 
evaluated using a mathematical model given by Adams (ref. 5). This value, along with the 
measured stick motion, provides the initial value for the deterministic pilot output u(t). 

The derivative of the pilot remnant is obtained by numerically differencing pilot-remnant 
values at consecutive points; the variable is evaluated from equations (6). The param- 

eter values are based on those given by Adams (ref. 5) and yield a stable system. For the 
variable variances, the same initial values used for case 1 of the simulated data were 

employed. Thus, the variances for the variables c(t), u(t), and ? 7 (t) were 0,00092 volt^, 

0.00046 volt^, and 0.00046 volt^, respectively; for the remaining variables, the value 0.1 was 
used. For the parameter variances, each standard deviation was set to 10 percent of the 
corresponding parameter value. The initial conditions for the first set of experimental data 
are given in column one of table IV. The estimates obtained after processing 1500 data 
points may be found as case 1 in the second column of table IV. The statistical parameters 
are in column one of table V. Most of the parameter estimates differ considerably from 
the corresponding initial values, particularly the negative values of Kj, and 

Further, these parameters represent an unstable system. The most revealing information is 
contained in the rms residual of the plant output, which is almost two orders of magnitude 
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larger than the known measurement noise standard deviation (0.0303 volt). The time histories 
of the stick-motion and plant-output residuals are given in figures 4 and 5. The behavior of 
the residuals in figure 5 suggests that the plant-output errors contain a signal as a result of 
errors in the mathematical model. Only the plant-output residuals are consistently larger than 
the measurement standard deviation whereas the corresponding variable c(t) standard devia- 
tion is small (0.0063 in table V). Therefore, it appears that the Kalman gain on this variable 
initially decreased rapidly as the measurement residuals decreased; the gain could not increase 
to compensate for later large measurement residuals. Since this phenomenon did not arise 
during the processing of the computer-generated data, the divergence of the plant-output 
residuals is a result either of processing very accurate measurements for a short period or of 
the presence of an unknown process noise (ref. 11). This situation is further aggravated by 
model biases caused by incorrect parameter values. The net result is an adjustment of the 
parameters; such adjustment continues until the system becomes unstable. 

As indicated by reference 11, all of the possible causes of the divergence may be 
treated by assuming the existence of a process noise of appropriate magnitude. The adaptive 
technique described in the appendix was chosen to calculate automatically the required process- 
noise statistics. Application of the original adaptive technique to this problem proved satis- 
factory. However, since the stick-motion residuals did not diverge in case 1, the modified 
technique described in the appendix was found to be sufficient to treat the plant-output 
divergence. In the modified technique only excessive plant-output residuals (that is, those 
exceeding the one-sigma value of 0.0303 volt) are used to adjust directly the Kalman filter 
gain for only the plant-output variable c(t). Changes in other gains and variables are the 
direct result of a normal application of the Kalman filter. In this way the modified tech- 
nique produces a compromise between the results of the Kalman filter and the original 

adaptive technique. 

The results of applying the modified technique are given as case 2 in tables IV and V. 
The rms plant-output residual is approximately the correct value while the rms stick-motion 
residual is little changed from that of case 1. Comparison of figures 4 and 5 with 

figures 6 and 7 further illustrates that only the plant output is significantly affected and 

that these residuals are more random. The new parameter estimates are clearly more reason- 
able than those in case 1. The system evaluated with these estimates is also stable. 

Table V indicates that by increasing the Kalman gain on only the plant output, many of 
the standard deviations have also increased to account for the uncertainty introduced into 
the modified estimation procedure by the fictitious process noise. 

Since the study using computer-generated noise indicated that iterating the Kalman filter 
over the data improved the parameter estimates, this iterative approach was also applied to 
the experimental data. The case 2 parameter estimates and variances are used to start the 
second iteration of the filter equations at 5.3 seconds. In the tables, case 3 gives the 
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estimates at the end of the third iteration. Although Ty 0, and appear 

to be converging to constant values, Tj appears to be decreasing rapidly. For this partic- 
ular problem all the parameters should be positive. The poor estimate of Tj undoubtedly 
results from the difficulty in separating the deterministic pilot-model parameters; a similar 
difficulty was observed in the case of the simulated data. Also, as seen previously, the 
observability of «2 highly questionable; in case 3, appears to be converging to a 

large positive number at best. The stability of case 3 is verified by the eigenvalues in 
table VI. 

Since the initial parameter values used in cases I to 3 should be good starting values 
for the estimation process, the only quantities which can be adjusted to improve the esti- 
mates of case 3 are the initial parameter variances. With this idea in mind, iteration of 
the Kalman filter was reinitiated using the original parameter values and new standard 
deviations equal to 3,2 percent of the parameter values. The results after three iterations 
are given as case 4 in tables IV and V. 

All of the parameters except have essentially converged to constants in case 4; 

a few additional iterations would only refine the values in the third or fourth place. F^or 
this case, all the parameters including Tj are reasonable; the smaller initial variances have 
reduced the tendency of any one parameter to diverge. Fwen appears that it would 

converge to a more reasonable, although large, value if enough iterations were attempted. 

As in the previous cases, the system evaluated with these parameter values is stable; the 
system eigenvalues are given in table VI. 

F.xamination of the correlations for case 4 shows that u(t) and pCt) have the 
expected large negative correlation (-0.97) and that only the correlations between u(t) 
and y^(t) (0.85), between Y^(0 and ? 7 (t) (-0.82), and between u(t) and (-0.81) 

are very significant. The (u(t),y^(t)j and |u(t),r 2 ^ correlatioiis are induced by equation (9); 
the |y^(t),77(t)j correlation is a result of the (u(t),77(t)) and |u(t),y^(t)j correlations. 

Comparison of the parameter standard deviations for cases 3 and 4 verifies the con- 
vergence of case 4. Since the initial standard deviations of cases 3 and 4 are, respectively, 

10 percent and 3.2 percent of the initial parameter values, the case 4 initial standard 
deviations are 32 percent of those of case 3. However, from table V, the parameter 
standard deviations of case 4 are 35 percent to 73 percent of those of case 3. The smaller 
standard deviations suggest that the standard deviations of case 4 are attaining steady values 
which are bounded below by the Cramer-Rao bound. This fact was further confirmed by 
examination of the behavior of the standard deviations on the third iteration; during this 
iteration the standard deviations changed much less than on previous iterations. 

The second set of experimental data was obtained with the same engineer as test 
subject; essentially the same closed-loop task was performed but with a different disturbing 
function. The measurement variances are identical to those for the first experiment; however, 
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since the disturbing function was different, the process-noise variance of the pilot remnant 
changed to 5.9 volts^/seconds . The initial conditions for the estimation problem, given 
in table VII, were obtained from the measured data in the same manner as they were com- 
puted for the first experiment. The state variances for each case are identical to the vari- 
ances chosen for the corresponding case of the first experiment. The initial parameters, 
which can be found in reference 5, yield a stable system. The initial parameter standard 
deviations are 10 percent of the parameter values. 

The results of case 1, in which the unmodified extended Kalman filter was applied 
to the data, are given in table VII and in figures 8 and 9. As with the previous data 
set, the plant-output residuals are excessively large and appear to contain a signal. This 
situation is apparently caused by the same errors that occurred in the case of the first 
experimental data set. Because of the erroneous parameter estimates, particularly the 
estimate of Kp the system has become unstable. 

Application of the complete adaptive technique to the data again indicated the use of 
the modified adaptive technique. Therefore, in case 2 one iteration of the modified tech- 
nique was used; this technique uses only the plant-output residuals and only the predicted 
variance of the plant output is directly altered. Figures 10 and 11 and the rms residuals 
in table VIII indicate that the divergence of the plant-output residuals has been corrected 
with little effect on the stick-motion residuals. Although the parameter estimates appear 
reasonable, the system of equations based on these estimates is unstable, and therefore 
inadequate to describe the physical system. 

As in the case of the previous data, the approach to preventing the system from 
becoming unstable is to restrict the changes in the parameter estimates by reducing the 
variances. In case 3, therefore, 3.2 percent of the parameter values was used for the corre- 
sponding initial standard deviations. The parameters given in table VII are the result of 
two iterations of the Kalman filter. Although these values appear reasonable, the system 
actually became unstable during the course of the second iteration. (See table VI.) Thus, 
the estimation procedure was abandoned for this set of standard deviations. 

In case 4, 1 percent of the parameter values was used for the initial standard deviations. 
The values in tables VII and VIII are the results obtained after four iterations of the Kalman 
filter. The system based on these values is stable. (See table VI.) With the exception 
of OL^ the estimates have essentially converged. Although the estimates and standard 
deviations of Tj and are sufficiently similar to suggest these parameters are highly 

correlated, such is not the case. The two parameters have very similar, but extremely small, 
correlations with several other variables and parameters; thus, since the correlation coeffi- 
cient for Tj and is only -0.0043, the two parameters only appear to be correlated. 

With the exception of the usual high correlation between u(t) and Tj(t), the greatest 
correlations are between u(t) and y^(t) (0.79) and between u(t) and (-0.86). 
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Therefore, the modified technique provides reasonable estimates of the state and of 
most of the parameters when applied to experimental data if a careful selection of the initial 
variances is made. In particular, reasonable estimates of the deterministic pilot-model param- 
eters can be made. 


CONCLUDING REMARKS 

An attempt has been made to estimate the parameters of a human pilot model by 
using the extended Kalman filter. The chosen human pilot model, originally written in 
transfer function form, was transformed to a differential-difference equation in the time 
domain; in this form the pilot model was amenable to treatment by the extended Kalman 
filter, although the time delay occurs implicitly. 

In studies utilizing computer-generated data, it was found that the deterministic pilot- 
model parameters, including the time delay, could be reasonably estimated if the magnitude 
of the process noise was sufficiently small, although the remnant to stick-motion ratio was 
realistically too small. The only difficulties appear to be separation of delayed term coeffi- 
cients into component parameters and correct estimation of a parameter in the pilot remnant. 

In the case of experimental data, it was found that the Kalman filter produced 
diverging plant-output residuals, apparently because of erroneous assumptions made on the 
mathematical model. Application of a modified adaptive technique was found to overcome 
this divergence problem. Careful selection of initial parameter variances produced reasonable 
estimates of the deterministic pilot-model parameters and maintained the stability of the 
system of equations. 

Therefore, the extended Kalman filter is recommended as a technique for estimating 
the parameters of the dynamical plant model and of the deterministic pilot model in a 
human tracking-task situation. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
July 9, 1975 


19 



APPENDIX 


AN ADAPTIVE TECHNIQUE THAT COMPENSATES FOR MODEL ERRORS 

In many estimation problems, errors in the state model arise from the use of simplified 
models or models containing erroneous parameter values. As a result of using such models, 
the estimation procedure fails to track the measurements properly and eventually diverges. 
Various techniques for alleviating this divergence phenomenon have been proposed (ref. 12). 
For the purposes of the current study a technique which could be easily automated was 
desired. 

The technique selected for this research is basically the scheme presented by Jazwinski 
in references 9 and 13. The philosophy of this technique is based on the fact that as long 
as the measurement residuals do not exceed the one standard deviation level, the Kalman 
filter is not altered. However, should the residuals exceed this level, a fictitious plant noise 
is computed from the measurement residuals. The fictitious noise is calculated in such a 
fashion that as the measurement residuals increase, the fictitious plant noise increases. The 
effect of the fictitious noise is to increase the predicted covariance matrix of equation (19) 
and, consequently, the filter gain. Thus, as the measurement residuals increase, the gain 
increases to “open” the filter. (That is, the filter uses more of the residuals to correct the 
predicted state.) If the residuals decrease to less than one standard deviation, the fictitious 
noise is not computed and the covariance matrix and gains reduce. 

With the use of this technique, it is necessary to rewrite two of the filter equations. 
Equation (18) is replaced by 


z(^k+lhk) y('-k|*-k) + f ^ ^ irzi^k^k)' Z (*-k-qhk-q)’ ^(^k-q)’ 

ytk 

+ r(k + 1) w(k + 1) (Al) 

where F(k +1) is a 13 X 1 coefficient matrix and w(k +1) is a zero mean 
Gaussian sequence with variance 
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The key to the adaptive process is to select the variance q(k + 1) (eq. (A2)) which 
reflects the previously stated philosophy. This selection is accomplished by estimating 
q(k + 1) from the measurement residuals; then q (k + 1) can be averaged over several 
measurement times to provide a statistically significant value. 


The variance q(k + 1) is estimated by maximizing the probability-density function 
of the residual with respect to q(k + 1); details of the derivation may be found in 
reference 9. The result of this derivation is that 


q(k + 1) = 


l_lj^2(k + 1). ejr^Ck + l)|q(k) = o| 


(If positive) 
(Otherwise) 


(A3) 



(A4) 


is tlic average iionnalized predicted residual at t^+i; 

^ Hj(«k+l) l’('ktll'k) "j^i'k+l) 




j=l 


[*^jj(^k+l) ^ *^ii(^k+l) 


1/2 


(A5) 


and 



(A6) 


In equations (A4) to (A5), mj|tj^^_]j and y(lk+lhk)j- indicate the jth compo- 
nent of the corresponding vectors, Rjj|t|,_^jj and the jth and ith diagonal 

elements of R^t|^+l), and Hj(tj^^jj, and Hj(tj^^lj are the Jth and ith rows of 
H(tk_|-i). In these equations, N is the number of measurements at time tj^-j-|; for the 
pilot-model problem, N = 2. The variance q(k + 1) was then averaged over the latest 
20 time points for this particular application. 

The major difficulty in applying this adaptive technique to a particular problem is in 
engineering an appropriate coefficient matrix F(k +1). Experimentation with the pilot 
model indicated the vector 


r(k + 1) = [0.707, 0, 0.5, 0, 0.5, 0, 0, . . ., O]'’^ 
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provided adequate results. However, since the results of case 1 based on the unmodified 
Kalman filter showed that only the plant-output residuals diverged, the use of the stick- 
motion residuals in the adaptive technique and the modification of the corresponding 
gains (i.e., the gains for the deterministic pilot output and the pilot remnant) should be 
unnecessary. Therefore, the technique was modified to use only the plant-output residuals; 
the coefficient matrix F(k + 1) was altered to influence only the plant-output estimate. 
The former modification requires restricting the summations in equations (A4) to (A6) to 
only those terms involving the plant output; the latter was accomplished by defining 
r(k -f- 1) by r(k -t- 1) = [0.707, 0, 0, 0, 0, . . ., 0]"^. These modifications were used 
in cases 2 to 4 for both sets of experimental data. 
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TABLE I.- INITIAL CONDITIONS FOR SIMULATED CASES 


Quantity 

True value 

Initial estimate 

(t = 0.25) 

(t = 0.25) 

C(t), V 

-0.9590 

-0.9496 

c(t), V/sec 

-5.6363 

-.5900 

u(t), V 

.6237 

0 

y4(t), V 

-.6880 

0 

T?(t), V 

.4105 

1.0109 

I 7 (t), v/sec 

2.5484 

2.48 

Kl 

1.0 

1.1 


.1 

.09 

r 

.1 

.12 

2 



7 

.5 

.6 

3 



Q 

.05 

.04 


4.427 

4.0 

“2 

10.0 

10.5 


50.0 

61.1 

^l"2 




-40.0 

-69.25 

^i\ ^ 1 / 
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TABLE II.- ESTIMATION RESULTS BASED ON SIMULATED DATA 



True 

Case 1 

Case 2 

Case 3 

Case 4 

Quantity 

value 

estimates 

estimates 

estimates 

estimates 

C(t), V 
c(t), V/sec 
u(t), V 
Y4(t). V 

T?(t), V 

V/sec 

(t = 10.25) 

-1.0692 

-3.0421 

2.2984 

-.9398 

.0734 

10.9367 

(t - 10.25) 

-1.0602 
-3.0225 
3.3659 
-2.4520 
-.9941 
1 1.8296 

(t = 10.25) 

-1.0602 

-3.0223 

3.2007 

-2.8760 

-.8289 

12.5505 

(t = 15.25) 

(t = 15.25) 

r 

1 

1.0 

1.1425 

1.1415 

1.0488 

0.9843 

.1 

.0865 

.0831 

.0936 

.0954 

.1 

.1 154 

.1060 

.1249 

.1435 

^ 1 
^3 

.5 

.6245 

.5782 

.5720 

.5282 

0 

.05 

.0399 

.0265 

.0401 

.0435 

“2 

4.427 

4.0014 

4.0650 

4.0980 

4.3534 

10.0 

10.4644 

7.6231 

10.4935 

10.5491 

K.r, 






1 ^ 

T r 
1 2 

50.0 

71.48 

74.93 

51.32 

37.98 

i 

1 






-40.0 

-82.15 

-81.84 

-57.27 

-46.81 



TABLE III.- STANDARD DEVIATIONS AND ROOT-MEAN-SQUARE 
(rms) MEASUREMENT RESIDUALS OF SIMULATED-DATA CASES 


Quantity 

Case 1 

Case 2 

Case 3 

Case 4 

Standard deviation 
of c(t) 

0.0063 

0.0063 

0.0063 

0.0063 

Standard deviation 
of c(t) 

.0201 

.0201 

.0201 

.0201 

Standard deviation 
of u(t) 

.3173 

2.5077 

.0580 

.1604 

Standard deviation 
of Y 4 (t) 

.5370 

4.9018 

.1584 

.4374 

Standard deviation 
of r?(t) 

.3187 

2.5079 

.0677 

.1628 

Standard deviation 
of T?(t) 

240.04 

240.01 

24.3396 

5.7711 

Standard deviation 
of Kj 

.0969 

.8868 

.0966 

.0648 

Standard deviation 
of Tj 

Standard deviation 
of r 

2 

.0081 

.0753 

.0082 

.0063 

.0106 

.1003 

.0108 

.0100 

Standard deviation 
of T 

3 

.0532 

.5038 

.0530 

.0368 

Standard deviation 
of d 

.0040 

.0353 

.0040 

.0056 

Standard deviation 
of 

.3999 

3.9432 

.3999 

.3983 

Standard deviation 
of 

1.0484 

10.2003 

1.0488 

1.0475 

rms residual, 
stick motion 

.00023 

.00002 

.00101 

.00654 

rms residual, 
plant output 

.0320 

.0320 

.0317 

.0240 



TABLE IV.- RESULTS BASED ON FIRST EXPERIMENTAL-DATA SET 


Quantity 

Initial 

values 

c(t), V ' 

(t = 5.3) 
-0.5 

c(t), V/sec 

-6.248 

u(t), V 

-.128 

y4(t), V 

.05 

T?(t), V 

.1 

f?(t), v/sec 

1.5615 

•^1 

.8 


.1 

1 


T 

.1 

2 


T 

.6 

3 


0 

.05 

"l 

4.427 

"2 

10.0 


Case 1 Case 2 
estimates estimates 
(t = 20.3) (t = 20.3) 

-9.9201 -4.0711 

-8.4231 -20.9579 

-.5683 -1.0824 

-.5711 -3.2801 

2.7453 3.2530 

5.5469 15.5910 

-.1463 .5661 

.3527 .0986 

.1669 .1034 

.0713 .4790 

.0544 .0706 

-.9218 3.4763 

-36.755 12.6328 


Case 3 
estimates 
(t = 20.3) 

Case 4 
estimates 
(t = 20.3) 

-4.0711 

-4.0709 

-21.0600 

-20.9791 

-1.3206 

-1.0170 

-5.3395 

-3.5209 

3.4880 

3.1870 

12.3238 

17.2188 

.4474 

.6354 

.0417 

.1099 

.1046 

.1094 

.4980 

.4855 

.0826 

.0655 

2.8634 

3.9206 

16.5106 

10.9260 
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TABLE V.- STANDARD DEVIATIONS AND ROOT-MEAN-SQUARE 
(rms) MEASUREMENT RESIDUALS FOR FIRST 
EXPERIMENTAL-DATA SET 


Quantity 

Case 1 

Case 2 

Case 3 

Case 4 

Standard deviation 
of c(t) 

. 

0.0063 

0.0297 

0.0297 

0.0297 ' 

Standard deviation 
of c(t) 

.0201 

.1002 

.0994 

.1006 

Standard deviation 
of u(t) 

.1638 

.2452 

.2006 

.1089 

Standard deviation 
of y^ft) 

.1243 

.4678 

.5326 

.2068 

Standard deviation 
of 77(t) 

.1664 

.2466 

20">2 

.1122 

Standard deviation 
of f?(t) 

5.2829 

5.4000 

5.3325 

5.2471 

Standard deviation 
of 

.0358 

.0682 

.0542 

.0226 

Standard deviation 
of r 

.0086 

.0092 

.0050 

.0030 

1 

Standard deviation 
of r 

.0025 

.0093 

.0086 

.0030 

Standard deviatic'ai 
of r 

3 

.0259 

.0528 

1 

.0488 

.0171 

Standard deviation 
of 0 

.0026 

.0016 

.001 1 

.0008 

Standard deviation 
of a j 

.3475 

.3290 

.2772 

.1235 

Standard deviation 
of 

! .9637 

.9197 

.8735 

.3063 

rms residual, 
stick motion 

.0116 

.0107 

.0101 

.0108 

rms residual, 
plant output 

2.8557 

.0339 

.0359 

.0340 
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I 


TABLE VI.- EIGENVALUES OF FIRST ORDER FADE MODEL FOR 
SELECTED EXPERIMENTAL-DATA CASES 


First data set. 

Fdrst data set. 

Second data set, 

Second data set, 

case 3 

case 4 

case 3 

case 4 

(-4.754, 0) 

(-29.363, 0) 

(-31.952, 0) 

(-36.575, 0) 

(-0.635. 2.624) 

(-15.780, 0) 

(-14.438, 0) 

(-15.615, 0) 

(-0.635, -2.624) 

(-3.191, 0) 

(-4.577, 0) 

(-4.159, 0) 

(-1.432, 3.803) 

(-0.229, 3.297) 

(0.026, 2.973) 

(-2.186, 2.310) 

(-1.432, -3.803) 

(-0.229, -3.297) 

(0.026, -2.973) 

(-2.186, -2.310) 

(-25.861, 6.886) 

(-1.960, 2.661) 

(-2.116, 2.546) 

(-0.130, 3.342) 

(-25.861. -6.886) 

(-1.960, -2.661) 

(-2.116, -2.546) 

(-0.130, -3.342) 
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TABLE VII.- RESULTS BASED ON SECOND EXPERIMENTAL-DATA SET 



• Initial 

Case 1 

Case 2 

Case 3 

Case 4 

Quantity 

values 

estimates 

estimates 

estimates 

estimates 


(t = 5.0) 

(t = 20.0) 

(t = 20.0) 

(t = 20.0) 

(t = 20.0) 

C(t), V 

-0.0155 

-2.4981 

-1.5373 

-1.5375 

-1.5376 

c(t), V/sec 

-4.6819 

-2.2670 

-4.2194 

-4.2693 

-4.3483 

u(t), V 

-.0811 

-.1618 

.0603 

-.0003 

-.0692 

Y4(t), V 

.3583 

.3185 

-.3228 

-.5054 

-.6790 

r?(t). V 

.05 

1.4936 

1.2697 

1.3307 

1.4012 

T?(t), V/sec 

.0 

3.8801 

5.6356 

5.9007 

5.9616 


.8 

-.1364 

.5488 

.6854 

.7629 

T 

1 

.1 

.1068 

.1132 

.1095 

.1034 

^2 

.1 

.0850 

.1147 

.1094 

.1035 

^3 

.4 

.2930 

.2944 

.3438 

.3818 

6 

.05 

.05568 

.0646 

.0613 

.0536 


4.427 

4.3639 

4.0958 

4.2329 

4.3721 


10.0 

.0199 

13.9415 

10.9624 

10.1136 


TABLE VIII.- STANDARD DEVIATIONS AND ROOT-MEAN-SQUARE 
(rms) MEASUREMENT RESIDUALS FOR SECOND 
EXPERIMENTAL-DATA SET 


Quantity 

Case 1 

Case 2 

Case 3 

Case 4 

Standard deviation 
of c(t) 

0.0063 

0.0246 

0.0246 

0.0247 

Standard deviation 
of c(t) 

.0201 

.0706 

.0714 

.0724 

Standard deviation 
of u(t) 

.0687 

.0331 

.0142 

.0053 

Standard deviation 
of Y 4 (t) 

.1626 

.0708 

.0344 

.0135 

Standard deviation 
of 7?(t) 

.0734 

.0418 

.0293 

.0262 

Standard deviation 
of r?(t) 

3.1752 

3.1645 

3.1599 

3.1580 

Standard deviation 
of Kj 

.0581 

.0716 

.0248 ' 

.0079 

Standard deviation 
of r 

1 

Standard deviation 
of T 

2 

.0093 

.0095 

.0031 . 

.0001 

.0096 

! 

.0095 

.0031 

.0001 

Standard deviation 
of T 

3 

.0373 

.0366 

! 

.0124 

.0040 

Standard deviation 
. of 0 

.0035 

.0024 

.0013 

.0004 

Standard deviation 
of Q!j 

.3357 

.3174 

.1332 

.0436 

Standard deviation 
of 0^2 

.9587 

.9337 

.3136 

.0997 

rms residual, 
stick motion 

.0164 

.0163 

.0165 

.0169 

rms residual, 
plant output 

1.3194 

.0373 

.0368 

.0360 





















Time, sec 


Figure 1 1 Plant-output error plotted against time for case 2 of 
second experimental-data set. 
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